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Abstract 

We propose a new solvable one-dimensional complex PT-symmetric potential as V{x) = 
ig sgn(x) |1 — exp(2|x|/o)| and study the spectrum oi H = —d?/dx‘^ + V{x). For smaller val¬ 
ues of a, <7 < 1, there is a finite number of real discrete eigenvalues. As a and g increase, there exist 
exceptional points (EPs), gn (for fixed values of a) causing a scarcity of real discrete eigenvalues, 
but there exists at least one. We also show these real discrete eigenvalues as poles of reflection 
coefficient. We find that the energy-eigenstates tpnix) satisfy (1): PTipnix) = l^hnix) and (2): 
FTtpEnix) = ipE^ix), for real and complex energy eigenvalues, respectively. 
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PT-symmetric quantum mechanics [1,22,24] which started theoretically [1-10,12] has pen¬ 
etrated well into the experimental and technological domains [11,14-16]. In this part of 
quantum mechanics, one considers non-real, non-Hermitian Hamiltonians which are invari¬ 
ant under the joint action of Parity (P : a; —)■ —x) and Time-reversal {T : i ^ —i). Even the 
most simple Hamiltonian H = — ^ + V{x) corresponding to the Schrodinger equation for 
these potentials has given astonishing results. Based on numerical computations. Bender 
and Boettcher [1] conjectured that the spectrum of VBB{x,e) = x^{ixY was entirely real 
when e > 0. This conjecture has later been proved [3]. Next, for —1 < e < 0 the spectrum 
consisted of a few real and the rest as complex-conjugate pairs of discrete eigenvalues. In 
the former case the energy eigenstates were also the eigenstates [1,2] of PT while the PT- 
symmetry was exact or unbroken. Interestingly, Vbb{x, 2) = — is a real Hermitian barrier 
(not a well), where the real positive discrete spectrum has been aptly interpreted [4] as the 
reflectivity zeros in scattering from flat-top potentials such as V{x) = — n = 1,2,3.... 

By complexifying Razavy’s real potential a quasi-exactly solvable potential was reported 
[5] displaying the phenomenon of broken and unbroken PT-symmetry. However, here P 
was taken as x ^ in/2 — x. A real Hermitian and a complex PT-symmetric Scarf H 
potential were found to have identical [6] spectra. Several other exactly solvable potentials 
were complexihed to produce [7] larger number of exactly solvable complex PT-symmetric 
potentials having real discrete spectrum. Existence of two branches of real discrete spectrum 
in complex PT-symmetric Scarf H was revealed [8] and interpreted in terms of quasi-parity 
[9]. Complex PT-symmetric Scarf H was shown to be an exactly solvable model displaying 
the spontaneous breaking of PT-symmetry when the strength of the imaginary part, IV 2 I, 
exceeded a critical value of 77-1-1/4, where —77(77 > 0) was the strength of the real part [10]. 
Such a phase transition of eigenvalues from real to complex conjugate pairs controlled by 
a critical parameter has inspired very interesting experiments in wave propagation [11] and 
optics where they realized PT-symmetry as equal gain and loss medium. Non-reciprocity of 
reflection [12] from such complex PT-symmetric mediums has given rise to novel phenomena 
like spectral singularity [13], coherent perfect absorption (CPA) without [14] and with [15,16] 
lasing, also see Ref. [17] for exactly solvable models of CPA of both kinds. 

Analyticity of a function in one dimension means continuity and differentiability in a 
domain. In a very interesting paper [18], the energy eigenspectrum of the potentials like 
V{x) = {ix)‘^\x\^ was studied to conjecture [18] that, except in rare cases (piecewise constant 


2 



potentials [19]), analyticity is an essential feature that is necessary for [a complex PT- 
symmetric] Hamiltonian to have the [entire] real spectrum. The words in square brackets 
were not mentioned earlier which we feel are necessary to be mentioned and emphasized. 
Following this, it was found that the real Hermitian potentials like V{x) = and |a;| 
having inhnite spectrum showed a scarcity [20] of real discrete eigenvalues when perturbed 
by V'{x) = ix\x\ and ix, respectively. It turned out [20] that if the total complex PT- 
symmetric potential is non-analytic, the real discrete eigenvalues may again be found but 
scarcely. 

In this Letter, we present a new exactly solvable complex PT-symmetric potential which 
is again non-analytic. We hnd both a hnite occurrence and a scarcity of real discrete eigen¬ 
values, depending upon the values of parameters. Scarcity in the spectrum is shown due to 
the presence of exceptional points (EPs). If 'ipn{x) is the energy-eigenstate, we hnd that (1); 
PT'^„(a;) = and (2) PT'ipE„{x) = [16,21] when the discrete eigenvalues are real 

and complex-conjugate pairs, respectively. The former (latter) is the property of un-broken 
(broken) PT-symmetry [16,21]. 

We consider the one-dimensional potential V{x) = ig sgn(a;) |1 — exp(2|a;|/a)| that can 
be rewritten as 

f ig[l — a; < 0 

V{x) = - (1) 

— 1], X > 0, 

without a loss of generality, let us assume that a,g > 0. When a tends to be large, V{x) —)■ 
2igx/a] this is so because ~ \ + z when \z\ is very small. The discrete spectrum of 
V{x) = i\x is known to be null [1]. 

For a; < 0, the time-independent Schrodinger equation can be written as 

^ + (p^ + = 0 , ( 2 ) 


where 


/2/i[E - ig] 

h 


1 -F i ^/2f^g 


This equation can be transformed into cylindrical Bessel equation [23], using y = 


Eq. (2) becomes 

y (iy2 + y^py~ ^^y + ^ ® ^ 

This Bessel equation [23] admits two sets of linearly independent solutions. For our purpose, 
we seek H^pl{y) and H^pl{y). Using the asymptotic property of cylindrical Hankel functions 
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[23], we find that 




~ \/2/7ry e 


[±i(j/-i^7r/2-7r/4)] ^ |^| ^ 


So the wave functions 


~ and ^ 2 ( 2 :) ~ 


—iRae, 


act as asymptotic forms of two linearly independent solutions of (1) on the left {x < 0). 
Seeing the dehnition of s in (3) we observe that ipiix) converges and 'ip 2 {x) diverges when 
X —)■ —00. Thus, for the bound state solution of (1), we write 

M^)= , ^< 0 . (7) 

Hipiisa) 

For a; > 0, by inserting the potential (1) in Schrodinger equation 

= 0 , ( 8 ) 


where q = Then using the transformation = sae^^'^, we get the modihed 

cylindrical Bessel Equation given below: 

- (z^ - ,/a^)iiz) = 0^ (9) 

LtZ (JjZ 

This second order differential equation has two linearly independent solutions I^{z),K^{z). 
Most importantly, we choose the latter noting that K^i^z) ~ when z ^ 00 . So, we 

write 

The solutions (7) and (10) are self-matched at a; = 0, by matching their derivative at a; = 0 
we get the energy eigenvalue formula 

f[E) = X%^(sa)K„a(sa) + i/“(so)A'',„(so) = 0, (11) 

for the real discrete spectrum of (1). For scattering states (positive energy continuum), we 
seek the solution of (4) when a; < 0 as 

ij{x) = ^ ( 12 ) 

ipi^x) = CKiqa{sae^^’^), a; > 0. (13) 
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By matching the solutions (12) and (13) at a; = 0, we obtain the reflection amplitude {B/A) = 

+ Kiqa{sa)H%^{sa) 


r{E) = i 


+ Kiga{sa)Hlp^ {sa) 


( 1 )' 


(14) 


The energy function, f{E), turns out to be the denominator of the reflection amplitude 
signifying that real energy poles of (14) are nothing but the real discrete eigenvalues of the 
exponential potential (1). The results (11) and (14) hold for g > 0. For < 0, in these 
results p and q; the Hankel functions and need to be interchanged. Eventually, 
eigenvalues of (1) are symmetric, in this regard En{—g) = En{g). 

The real discrete spectra of (1) can be obtained by the real zeros of f{E) = 0, however, 
finding them by locating the real discrete energies, where the reflection probability {R{E) = 
|r(E)p) becomes infinity, is more interesting and presentable (see Figs. 1). Usually, the 
real energy, where reflection, R{E), and transmission, T{E), probabilities become infinite is 
called a spectral singularity [13] for a scattering potential which is such that U(±cx)) = 0. 
Here, as the potential diverges to ±fcx), the spectral singularities are the positive energy 
bound states. A recent study [22] of the common negative-energy singularities of R{E) 
and T{E) is shown to yield the real discrete spectrum of three parametric domains of the 
complex Scarf II which is neither FT symmetric nor (apparently) pseudo-Hermitian. 

In every numerical calculation in the sequel we shall be using 2p = 1 = . In Fig. 1, we 

plot R{E) foTg = l and a = 0.5,1,2,5 in (a,b,c,d), respectively. The poles in R{E) at dis¬ 
crete real energies (eigenvalues: roots of Eq. (11)) appear as sharp peaks, just after the peaks 
the minima indicate coalescing of a pair of real eigenvalues into complex- conjugate pairs. 
The number of eigenvalues reduces from 9 to just 1 as we increase a from 0.5 to 5 in Fig. 1. In 
Fig. 1(b), we see five clear spikes at real energies E = 3.27651(Eo), 8.83705(Ei), 13.7572(772), 
21.3301(773), 25.6883(774) and one wide maximum at 77 = 37.5832 ± 2.6879f(775 e). We And 
that the bound eigenstates for real energies are like 'ip{x) = (pei^x) + ipoix), (f>e,o{x) are even 
and odd parity functions satisfying (;/)e,o(±C)o) = 0 [24]. In Fig. 2, we plot real and imaginary 
parts of 'ip{x) for E = Eq and 77i. The even (solid line) and odd (dashed line) parities of 
these parts testify to 

F^^/Jnix) = = +lpJn{x), (15) 


signifying that PT-symmetry is exact (unbroken). The next E^ and Eq the poles of 77(77) 
or the solution of /(77) = 0 are the complex-conjugate pairs of energy eigenvalues. The 
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interesting signature of the spontaneous breaking of PT-symmetry results in the loss of 
dehnite parities of the real and imaginary parts of the eigenstates PT('^(a;)) or ip{x) and 
both differ from each other (see Fig. 3). In the case of spontaneous breaking of PT-symmetry, 
the flipping [16,21] of eigenstates takes place as 

PT^jJE{x) = PT'ipE* = aeTZ, (16) 

For the present model, see Fig. 3, for the eigenstates 'ip 5 {x),'ipQ{x) corresponding to complex- 
conjugate eigenvalues. Figures (a) and (b) are mirror images of each other about y—axis and 
hgures (c) and (d) are minus times the mirror image of each other about y—axis. Thus, in 
the case of the exponential potential (1), we have a = 0 in (16). 

For a = 1, we hnd three consecutive exceptional points (values of g) as 0.74, 1.58, 5.35. 
The least number of real discrete eigenvalues is 7 when g < 0.74. Exact numbers of real 
discrete eigenvalues for 0.74 < g < 1.58, 1.58 < g < 5.35 and g > 5.34 are 5, 3 and 
1, respectively (see Fig. 4). Additionally, for g = Q.l and g = 0.2, we get 12 and 9 
eigenvalues, respectively. Similarly, for a = 0.5, we hnd four consecutive exceptional points 
as g = 1.62,2.96,6.29,21.38 (see Fig. 5). When g < 1.62 there are at least 9 eigenvalues. 
Then for 1.62 < g < 2.96, 2.96 < g < 6.29, 6.28 < g < 31.38, in Fig. 5, one can see that 
there are 7, 5, 3 and 1 eigenvalues, respectively. When the parameters a and g are very small 
there may be an abundance of real eigenvalues but the entire discrete spectrum cannot be 
real. We conjecture that in the proposed non-analytic complex PT-symmetric exponential 
potential the entire discrete spectrum cannot be real. For higher values of a and g, there will 
be scarcity of real discrete eigenvalues; nevertheless there will be at least one real discrete 
eigenvalue but for the limiting case of a —)■ cx) (when V{x) (1)) it passes over to i\x. 

The behaviour En{g) for g G {g 2 ,g 3 ) in Fig. 4 and for g G {g 3 ,gA) in Fig. 5 where there 
are only three real discrete eigenvalues, is akin to that of the quasi-exactly solvable complex 
PT-symmetric potential for M = 3 [5]. There the ground state is un-paired with Eq = 5 — 
but Ef = 7-z‘^± 2Vl - 4;^^ 0<z< 1/2. 

Strangely, despite having analytic eigenstate given by Eqs. (7) and (10), it is indeed 
challenging to prove elegant properties (15) and (16) of energy eigenstates in PT-symmetry. 
Just to bring out the level of difficulty, we at best can express [23] 

Hl^\z) = + 1/2, 2z/ + 1; -2iz), 

V TT 

K^{z) = V2e-^{2zYU{iJ + l/2,2iJ + l-,2z). (17) 
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Here U{a, b] z) is second Kummer’s function which is solution of confluent hyper-geometric 
equation with limited properties [23]. The argument of these functions in our case is like = 
(si-|-is 2 )e^/“ under PT-symmetry it will change to z = (si —Such a transformation 
of U does not seem to be available. 

Recently, it has been shown that square integrable eigenstates of a complex potential 
satisfying Dirichlet boundary condition ■^(icx)) = 0, regardless of real or non-real discrete 
eigenvalues, are orthogonal as [22] 

/ OO 

^/Jmix)^/Jnix)dx = 0. (18) 

•OO 

We hud that the same holds true here as well. For instance, 6 eigenstates corresponding to 
the case when g = I and a = 1 (Figs. 1(b), 2, 3) are mutually orthogonal as (18). In case 
we use PT-scalar product as < 'iprn{—x)\'ipn{x) > [2,9,10], then due to the property (16), 
'ip 5 {x) and 'ipei^x) corresponding to complex-conjugate eigenvalues are self- orthogonal and 
they have PT-norm as zero. 

We conclude by asserting that we have presented a new solvable complex PT-symmetric 
non-analytic exponential potential which has both a finite occurrence and a scarcity 
(existence of exceptional points) of real discrete spectrum depending on the values of 
its parameters. Based on our calculations, we conjecture that this potential cannot 
have entire discrete spectrum as real. In other words, there is no parametric separation 
of unbroken and broken PT-symmetry. There exists exceptional points (EPs) (several 
critical values of the parameter g for fixed values of a) below which there are only a hnite 
number of real discrete eigenvalues. However, the states having real eigenvalue are also 
degenerate eigenstates of PT all having eigenvalue of PT as 1. It may be recalled that 
for PT-symmetric Scarf H, these are (—1)”. Next, in the case of complex-conjugate pairs 
of eigenvalues, the action of PT flips the eigenstate 'ipE^x) to %Ije*{x). This potential has 
at least one real eigenvalue for hnite values of a. An interesting anti-climax of the new 
solvable model is that the commonly known properties of Hankel and modihed Bessel 
functions fall short to prove the elegant properties (15,16) of eigenstates under the action 
of PT; we have conhrmed them numerically. This physical requirement may inspire one to 
look for some new or rarely known properties of these aforementioned higher order functions. 
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FIG. 1: Spikes in reflectivity, R{E) (14) indicating real discrete eigenvalues (solutions of Eq.(ll)) 
of the newly proposed potential (1). Here g = 1 and a is being varied: (a) 0.5, (b) 1.0, (c) 2.0, 
(d) 5.0. In (a,b,c), see a hump after discrete energy poles indicating the next two eigenvalues have 
merged to become complex-conjugate pairs. One can actually choose a value of g for which the last 
two discrete energies will be very close by as real doublets. Notice that as a increases, the number 
of real eigenvalues reduces from 9 to just 1. In fact, for very large values of a there will be no real 
eigenvalue. This is so because the potential (1) passes over to V{x) = iXx which is known to have 
discrete spectrum as null. 


10 















(a) 

^0 


(a) 




FIG. 2: The typical complex eigenstates corresponding to real eigenvalues (a) Eq = 3.27651 and 
(b) El = 8.83705. The solid (dashed) lines show the real (imaginary) parts. Here a = 1 and 
S' = 1 as in Fig. 1(b). In general, we find that these eigenstates are square integrable of the type: 
= 4'e,nix) + i(j)o,nix), denote even(odd). 


(a) 




(b) 



'-A 1 


(d) 



FIG. 3: Typical eigenstates y^{x) and y%{x) corresponding to discrete complex-conjugate energy 
eigenvalues, 37.5832 ± 25.6883z, respectively. Potential parameters are same as that of Figs. 1(b) 
and 2. Notice that these eigenvalues appear in Fig. 1(b) as a hump after 4 real discrete energy peaks 
(poles). One can closely look into these figures to make out the acclaimed result, i.e., PT(V^ 5 (x)) = 
^^{x). More directly, see that the figures (a) and (c) are reflection of each other about y—axis, 
but (b) and (d) are minus times the mirror reflection of the each other. Numerically, '^ 5 ( 1 ) = 
0.661638 0.121078i,V'5(-l) = 1.02862 - 0.201041i and V’6(l) = 1.02862 + 0.201041i,'06(-l) = 

0.661638 - 0.121078i. 
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FIG. 4: Real discrete energy eigenvalues of the exponential potential (1) as a function of g for 
a = 1.0; there are three EPs: gi = Q.7A,g2 = 1.58 ,173 = 5.35. 




FIG. 5: Real discrete energy eigenvalues of the exponential potential (1) as a function of g for 
a = 0.5: there are three EPs: gi = 1 . 62 , 5:2 = 2 . 96,53 = 6 . 29,54 = 21.38. 
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